Statistics of speckle patterns 
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We develop a general method for calculating statistical properties of the speckle pattern of coherent 
waves propagating in disordered media. In some aspects this method is similar to the Boltzmann- 
Langevin approach for the calculation of classical fluctuations. We apply the method to the case 
where the incident wave experiences many small angle scattering events during propagation, but 
the total angle change remains small. In many aspects our results for this case are different from 
results previously known in the literature. The correlation function of the wave intensity at two 
points separated by a distance r, has a long range character. It decays as a power of r and changes 
sign. We also consider sensitivities of the speckles to changes of external parameters, such as the 
wave frequency and the incidence angle. 

PACS numbers: 72.15.Rn, 73.20.Fz, 73.23.-b 
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In this article we consider statistical properties of 
waves propagating trough a disordered medium and de- 
scribed by the stationary (scalar) wave equation, 



(v2 + fc2ri2(r)) =0, 



(1) 



where n(r) is the index of refraction assumed to be a ran- 
dom Gausian function. This problem is relevant for a va- 
riety of important physical situations, ranging from elec- 
tromagnetic waves propagating through the interstellar 
space or the atmosphere, to electron transport in disor- 
dered conductors. The wave density, /(r) = |'(/'(r)p, ex- 
hibits sample specific random fluctuations (speckles) due 
to interference of waves traveling along different paths. 
The statistical properties of speckles have been discussed 
in the past. The problem can be characterized by sev- 
eral characteristic lengths: the propagation distance, Z , 
the elastic mean free path, H. and the transport length 
(■tr ~ ^7^*0, which is the typical distance for backscat- 
tering. Here is the typical scattering angle on 

the distance £, and ^ is the correlation length of n(r). 
In the diffusive regime Z ^ £tr the problem has been 
studied in Refs. The hmit of "directed waves" 

itr ^ Z £ has been studied in many papers, see for 
example 0,0,00 and references therein. In the latter 
case the wave experiences many small angle scattering 
events, but the total change of its propagation angle 9 
remains small. 

In this Letter we develop a general method for calcu- 
lating speckle correlations, which enables us to treat both 
diffusive and directed wave cases on equal footing. It is 
similar, in some of its aspects, to the Langevin scheme 
for the description of classical fluctuations 8, (oj, 0| . We 
shall demonstrate the method for the case itr ^ Z ^ £ 
by calculating the speckle correlations and their sensi- 
tivity to various perturbations, such as a change in the 
frequency of the wave, its incidence angle and a change 
of the refraction coefficient. In many aspects, our re- 
sults differ from those obtained in the previous stud- 
ies 0, IE 0) 0] • Among the differences are the slow power 
law decay of the density correlator as a function of coor- 




FIG. 1: The asymptotic behavior of the density correlation 
function C(p) — {SI{p)5I{0)). X is the light wavelength, £ is 
the elastic mean free path, Z is the slab width, and 9o and 9 
are the typical scattering angle of a ray traveling a distance £ 
and Z respectively. 



dinates and its change of sign, see Fig. ^ 

The central object of our approach is the ray distribu- 
tion function /(r,s) 



/(r,s) = 



r 

^+2 



(2) 



which is the probability of finding a ray at point r point- 
ing in the direction specified by the unit vector s, in par- 
ticular J(r) — J c?s/(r,s). The average distribution func- 
tion (/(r, s)) satisfies the Boltzmann kinetic equation. 



s-^^ = /.{(/(r,s))}. 



(3) 



dr 

IstUf)} = I <f~sW{s-s) [(/(r,s)) - (/(r,s))] .(4) 



Here (...) denotes averaging over the random realizations 
of n(r), the integral over the directions is normalized to 
unity, / (Ps = 1, and W{5s) = ^ / (Prg{v)e'^^^ '' is the 
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FIG. 2: Diagrams describing wave density correlations. The 
solid lines denote the disorder-averaged Green functions of 
Eq. Q and the dashed lines denote the disorder correlator, 
4k'^{5n{r)Sn{r')). 



probability, per unit length of propagation, for changing 
the ray direction by ds, and g{r — r') = (n(r)n(r')) — 
(n(r))(n(r')) is the disorder correlation function. 

Correlations of the fluctuations, Sf = f — (f), may be 
evaluated using the Langevin-type equation, 

s-^~Ist{Sf} = C{r,s), (5) 

where £ is a random Langevin source, with a zero mean 
and the variance given by. 



(/:(r,s)£(r',s')) = ^(5(r-r') 



<5(s-s')(/(r,s)) 



d2siW^(s-Si)(/(r,Si)) - (/(r,s))W^(s-s')(/(r,s')) 



(6) 



r 



One can prove Eas. (1316(1 using the standard impurity 
diagram technique 11]. Equations H3I4|I are obtained 
by summing the ladder diagrams in Fig. El a), whereas 
Eqs. (|5I()|) follow from diagrams b)-d) in Fig.|21 

Equations I|3I6() are valid when the mean free path 
is sufficiently large, £ = {J ds'^W{s))~^ > £,'^/X, , and 
|r — r'l ^ A. Here r and r' are observation points, 
? — [J (Prr'^ g{r)/3 J d'^rg{r)]^/^ > A is the disorder cor- 
relation length, and A — 27r/fc is the wave length. For 
|r — r'l < A Eqs. H3I6|) are not valid and to evaluate 
the correlation function one has to calculate the diagram 
shown in Fig. El e). 

In the regime of small angle scattering, ^ 1: and 
when |r — r'l 3> ? the rays undergo diffusion in the space 
of directions, s. In this case Eqs. (|3m|l reduce to a set of 
angular diffusion equations, 



g(/(r,s)) 

dr 

dSf{r,s) 

® a 

or 

(jf(r,s)j^(f,s)) 



DeV2(/(r,s)>, 



(7) 



V, [i?eV,,5/(r,s)-ji(r,s)], (8) 
2nDe{f)^ 



fc2 



-6^pSis~s)Sir-r). (9) 



Here j^(r, s) are the Langevin current sources, Dg = 
i^-i = i / d^s'{l - s • s')W^(s - s') is the diffusion con- 
stant in the space of angles s, and = 9-^ + jj^^^ 

is the gradient operator, with (f> = (— sin0, cos 0, 0), and 

6 — (cos (f) cos 9, — sin (j) cos 6, — sin (/>) . 

Further simplification emerges at larger spatial scales, 
|r — r'l itr- In this case Eqs. H3l6f) can be reduced to 
diffusion equations V^{I) = 0, and V{-DV5I + 

J) = 0. Here D = Itr/^ is the (real space) diffusion 
constant, and correlation function for Langevin currents 
has the form: {J^{v)Jp{T')) = ^{I{r)yS^f,S{r-r'). In 
this case, the correlation function {6I{r)SI{r')) decays as 
l/r and as 1 /r^ for r ^ itr and X <^ r ^ £tr respectively. 

The scheme presented above can be generalized to 
treat speckle sensitivities to changes of external param- 
eters, such as the wavenumber, Afc, or smooth changes 
of the refractive index, An(r). We characterize the sen- 
sitivities by the correlation function ((5/(0)5/(7)), where 

7 — Afc -I- fcAn(r). In this case the Langevin source cor- 
relator is given by. 



(£(r,s;0)/:(r',s';7)) = -^<5(r- 



0^ ,5(s-s')/,(r,s) / d2siM^(s-Si)/_,(r,Si)-/,(r,s)I^(s-s')/_,(r,s') 



(10) 



where /±(r, s) satisfies the equation, The method based on Eqs. is similar to the 

Langevin approach describing classical time and space 

df±(r,s) fiuctuations of a single particle distribution function 

s ^- /st{/±(r,s)} =±i7/±(r,s). (11) 
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[M ^1 ■ The fundamental difference between our prob- 
lem and that of classical fluctuations manifests itself in 
the form of the correlation function of the Langevin 
sources. In the classical problem the Langevin sources 
are ^-correlated both in time and space and their vari- 
ance is proportional to the average distribution function 
(/). In contrast, in Eqs. H6I10() the Langevin source vari- 
ance is quadratic in (/) and (5-correlated only in space. 

To illustrate the use of Eqs. H3I1HI we consider the case 
when a wave of intensity Iq is incident on a disordered 
slab of thickness Z, such that £ <^ Z <^ £tr, as shown 
in the inset of Fig. ^ The results presented below are 
calculated to leading order in the small total scattering 
angle, 6^ = DgZ. In this case the correlation function 
C(r — r') = ((5/(r)(5/(r')) is strongly anisotropic (here 
61 (r) = /(r) — (/(r))). Therefore below we shall use the 
notation: r = {z,p), where p denotes a two-component 
vector in the plane perpendicular to the z-axis and z de- 
notes the distance between observation points along the 
z axis. When p = 0, i.e the observation points are located 
along the z axis, the correlation function for z ^ Z is, 



C{z) 



^0 



(12) 



4fc26»4z2 ■ 

Equation 112|l matches the results for the diffusive case 
J»i I Z ^ Itr when 9 is of order unity. 

When z < p/9, i.e. the observation points are located 
essentially on a plane perpendicular to the z axis, a gen- 
eral formula for C{p), can be derived from Eqs. (|7I9|I . 



Cip) 



-'o 



d 

X — exp 
dC, 



C-z 

2 rC 

' £ 



dqqJoiqp) 



, (13) 



where g{p) ~ J dzg{\/ p^ + z^) / J dzg{z). The integral 
in Eq. H13|l contains a term proportional to a (5-function, 

2Dek^z ^(p)- This term represents the rapidly decaying 
(at p ~ A/6') part of the correlator and corresponds to the 
part of diagram Fig. |21 b) without the impurity ladders 
after the Hikami box, see diagram e). The (5- function 
term results from the semiclassical approximation em- 
ployed in the derivation of Eqs. H3I11(I . which limits the 
spacial resolution to (5p 3> A. In order to resolve the spa- 
cial structure of the short distance part of the correlator 
diagram e) needs to be evaluated more accurately. This 
gives the following asymptotic behavior, 



C(P) 



^-2{kBpf 
bi 



if 
if 

if 

if 



pr^a\/6, 

^6io < p < ez, 

0Z«p« ^, 



(14) 



where — log{k£6^ / 9q) , bi = dxg{x) is a con- 
stant of order unity, 62 = 3^/^r(5/3)/8 0.163, and 



63 = 27/128 « 0.21. The tail of the correlation function 
(the regime p > Z9'^/9o) is also described by Eq. ifT^ 
and depends on the precise form of the disorder correla- 
tion g{r), since this limit is dominated by rare scattering 
events. The qualitative form of the function C(/9) is shown 
in Fig.HI 

Let us consider now the statistics of density, inte- 
grated over a disk of radius R, P = Jp^ji d^pl{p)- Using 
Eqs. (O in we get. 



{{5Pf 



i?4 



n 

2WW 
h' -5- 



k'^e^Ba ' 



^6io<i?<6'Z, (15) 



where b'^ = 26i7r/3, b'^ = 3'*/3r(5/6)7r/2"/3r(7/6), and 
6(, = y372^. 

Consider now the sensitivity of the integrated density 
P{uj) to a change in the wave frequency Aw = cAfc, 
where c is the speed of the wave. It can be character- 
ized by the experimentally accessible quantity. 



((P(c^ + Ac.)-P(c.))^) 
{{5Pf) 



Alo 



where 



UJ 



2 e^z' 



-r ' (16) 



(17) 



A qualitative explanation of the scale w* is similar to that 
given for the sensitivity of the conductance fluctuations 
[1^ Hsf . Let us estimate the characteristic change in 
the phase of a typical orbit due to the frequency change 
Auo : The typical length spread of the orbits is of or- 
der O^Z. Therefore the phase difference is AkZO^ where 
Afc = Aw/c is the change in the wavenumbcr. Thus a 
complete change of the speckle pattern occurs when the 
phase, AioZfP jc is of order one, namely Auj ^ c/ZO^ , in 
agreement with Ref. Q. 

As another application of our scheme let us con- 
sider the sensitivity of speckles to a change in the in- 
cidence angle, </>, of the wave, see the inset in Fig. 2] 
(thus the incident wave function now has a form -0 — 
-y/To exp [ifcz cos 0-1-1 fcp sin (/)]). One may characterize this 
sensitivity by the correlation function. 



{{SPf) 



1 

'8 



k^£'^Dl 



(18) 



where 0* — {6kZ)^^, /3 is a factor of order unity, and it 
is assumed that R <C OZ. The first term of this equation 
follows from formula Q and the correlator IjlOl) . where 
f± satisfies equation Hll(l with initial conditions 



f±{p,z = 0) = loe 



S{s - So), 



(19) 
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where Sq = {cos(j),s±) « (l,s^), with |sj_| = sincj) « (j), 
assuming 0^1. The second term in the right hand side 
of Eq. ifT^ is computed from diagrams of the type shown 
in Fig. 121 d), containing two Hikami boxes. It represents 
a small correction in the parameter S,/£9o, however, this 
term becomes the dominant contribution when 3> 0* . 

On a more general level Eqs. (|3I11() provide a frame- 
work for calculation of time correlation of speckle pat- 
terns in the case where the refraction index changes in 
time, provided these changes are much slower than the 
time of propagation of the wave to the observation point. 

Our results also may be easily extended to cases with 
light polarization, optically active media, Faraday effect, 
and coherent short wave pulses as long as their duration 
is longer than r = £/c. These issues are left for future 
studies. 

The problem considered here is similar to the problem 
of universal conductance fluctuations in metallic samples 
[l^ which are also of interference nature. There- 
fore we would like to discuss the relation between the 
two problems. In the single particle approximation the 
conductance of a metallic sample, G ~ / d(f>T{(f>) can be 
expressed in terms of the electron transmission probabil- 
ity through the sample, T{(f>) cx J dpd'^s (z • s)/(p, s), in- 
tegrated over the incidence angle of the incoming wave, 
(j). Thus the variance of the conductance fluctuations, 
6G is proportional to a double integral of the correlation 
function {6T{(t))6T{(p')) . In principle, the latter can be 
calculated using Eqs. H3I11() . or, cquivalently, by calcu- 
lating diagrams shown in Fig. El b)-d). However, as we 
explain below, this does not account for the conductance 
fluctuations. The latter arise from diagrams of the from 
shown in Fig.|21f). 

In the limit of directed waves, Z <g; itr-, there is no 
backscattering. Therefore the transmission probability 
does not fluctuate, ST{(j)) ^ SG = 0. Thus to com- 
pare the two problems we have to consider the diffu- 
sive case, Z ^ Itr, where 6G ^ e^/h. The correlation 
function {dT{(j))ST{ip')) , in the diffusive regime, still has 
a structure similar to the correlation function given by 
Eq. H18() 3J. Namely, it contains two contributions. The 
first contribution comes from diagrams (21 b)-d) , and de- 
scribes relatively strong fluctuations of the transmission 
coefficient. However, it is very sensitive to the change of 
(j), and after the integration over gives a small contri- 
bution to SG^. The second contribution originates from 
diagrams of the type Fig. Of), and is analogous to the sec- 
ond term in Eq. (|18|l . Although its amplitude is smaller 
than that of the first term, it is insensitive to the change 
of (f>, and after the integration over yields the dominant 
contribution to SG^. Thus, conductance fluctuations are 
not described by Eqs. (|3I11|I . and should be calculated 



from the diagrams of the type shown in Fig. [21 f) (see the 
corresponding discussion in Ref. [3j). 

Finally we would like to mention that the results pre- 
sented above substantially differ from those known in the 
literature (see for example Refs. 0,0, IHH)- First, the 
correlation function H13() exhibits a universal long range 
power law behavior in a wide range of values of p. The 
only non-universal regimes are at the tail, p S> Zd^ /Oq 
, and the short distance region, p In contrast, in 

the results presented in Refs. 0,|EIE0 C{p) depends on 
the detailed form of g{r), and usually decays exponen- 
tially at p > ^. Second, in contrast with previous results, 
C{p) changes its sign as a function of p, which is a con- 
sequence of the current conservation. This conservation 
law also implies, that the fluctuations of the integrated 
intensity over disks of radius R > Z9 is proportional to 
R, see Eg. IIISI). rather than R^, as would follow from 

Refs. miai. 
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